Fate of chiral critical point under the strong isospin asymmetry 



Department of Physics, 



o 

(N 

C3 



43 

Oh- 
i . 

Oh- 

D ■ 

43" 



> 

O 
On 

o 



x 

S3 



Hiroaki AbukG 

Tokyo University of Science, 

(Dated: May 2, 2013) 



Tokyo 162-8601, Japan 



We study the influence of the isospin asymmetry on the phase structure of strongly interacting 
quark matter near the critical point (CP) using a Ginzburg-Landau approach. The effect is found 
to be drastic, not only bringing about the shift of the location of the CP, but resulting in a rich 
phase structure in the vicinity of the CP. In particular, new tricritical and triple points emerge 
as soon as the isospin density becomes finite. Moreover, we find the CP being washed out from 
the phase diagram due to the stabilization of a homogeneous charged pion condensate when the 
isospin chemical potential exceeds a critical value. We derive a model-independent universal relation 
between the critical isospin chemical potential and the chiral condensate at the CP. We also study 
the effect of the U(1)a anomaly on the phase transition to the pion condensate in the vicinity of 
chiral crossover. 



PACS numbers: 12.38.Mh, 21.65.Qr 



I. INTRODUCTION 

The phase diagram of QCD at finite temperature 
and/or finite density is the subject of extensive theo- 
retical and experimental studies. In particular, several 
approaches to QCD with two light flavors suggest the 
existence of a critical point (CP) at which the first-order 
chiral phase transition turns into a crossover [l| . Despite 
many efforts based on the first principle calculations [2j, 
not only the precise location of the CP, but its existence 
itself remains controversial. 

In the chiral limit with vanishing quark mass, the CP 
becomes a tricritical point (TCP). The effect of a finite 
quark mass is, thus, rather simple; just to smear the 
second-order chiral phase transition to a crossover and, 
accordingly turn the TCP into the CP. Our focus here is 
the other important ingredient in realistic systems, the 
effect of an isospin asymmetry. Such a flavor symmetry 
breaking is caused by a neutrality constraint that should 
be imposed in any bulk systems to prevent the diverging 
energy density. 

The isospin imbalance is known to bring a rich vari- 
ety of color superconducting phases at high density [3|. 
On the other hand, QCD at large isospin density was 
first studied in Ref. J4J, and it was shown that the QCD 
vacuum develops a pion condensate (PIC) as soon as 
| /ii| > m^, where m w and /ii are the vacuum pion mass 
and the isospin chemical potential, respectively. The PIC 
can be viewed as a relativistic superfluid that exhibits a 
crossover from a Bose-Einstein condensate (BEC) of pi- 
ons to a superfluidity of the Bardecn-Cooper-Schrieffer 
(BCS) type [5]. So far, several model analyses have been 
made for the PIC at finite temperature and/or quark 
density |5l-tl3|. However, to our knowledge, there is at 
present, no systematic analysis based on the Ginzburg- 
Landau (GL) approach focusing on the isospin effects on 
the CP. This is what we present here for the first time. 



Our GL framework is advantageous to other ap- 
proaches in the sense that it can give model-independent 
predictions near the CP. Since we are interested in the re- 
sponse of the CP and phases in its neighborhood against 
nonzero /ii, our strategy is to take fj,i as a perturbative 
field and expand the GL functional with respect to it. We 
use a quark loop approximation to reduce the number of 
GL couplings. This approximation should be valid, in 
particular, if it is located at a large fugacity region. 

In this paper, we restrict the analysis to homogeneous 
phases only, leaving more detailed analysis to future work 
14 1 . This is, to some extent, an extension of our previous 
work [15| to the situation off the chiral limit introducing a 
finite quark mass. This paper is organized as follows. In 
Sec.|TTl we derive a general GL potential up to the fourth 
order in fields and discuss the effects of fj,\ and U(1)a 
breaking at the vicinity of chiral crossover. In Sec. lIIIl we 
extend the GL potential up to the sixth order to discuss 
the isospin effect on the CP. Based on this, we clarify how 
the CP and its neighborhood are affected by the inclusion 
of isospin asymmetry. In Sec. IIV1 we summarize. 



II. GINZBURG-LANDAU APPROACH AT 
FOURTH ORDER 

Let us start with writing the most general GL potential 
for the two chiral four- vectors cf> = (a, tz) and its parity 
partner (p — (r]',a) with a ~ (<jq), t ~ (qi"f 5 Tq), r( ~ 
(qijsq), and a ~ (q~Tq). At the fourth order in <f> and ip, 
the chiral [SU(2) L x SU(2) R ~ 0(4)] symmetric part of 
the GL potential should take the form 
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where (</>, ip) = arf+a-ir is the inner product. If a 2 = a' 2l 
the potential possesses the additional U(1)a symmetry, 
which is actually violated in QCD via the axial anomaly, 
so, typically, a 2 ^ a' 2 - The current quark mass and the 



isospin chemical potential add 0(4) noninvariant terms 
to the potential. At the second order in fields, 



da, 



SB 



■ha + Pi(<Ta 3 + T]'m)+02ir±, (1) 



where ir±_ = (^1,^2) is the charged pion doublet. The 
first term is due to the current quark mass, and it ex- 
plicitly breaks 0(4) down to SU(2) V ~ 0(3). The sec- 
ond and third terms are due to the finite isospin density, 
which violate the isospin 0(3) into Ui 3 (l) ~ 0(2), the 
rotation about the third axis of isospin space. The GL 
coupling h is proportional to the quark mass m for light 
flavors. On the other hand, since the operator in the 
second (third) term is even (odd) under the isospin flip 
u «-» d, we have 0\ oc /ii and 0% ex /i 2 at the leading order 
in expansion in \x\. When \i\ ^ 0, and h 7^ 0, we may 
anticipate the realization of following two phases: 

(i) The chiral symmetry broken phase (xSB): the 

phase with a 7^ 0, which might be accompanied by a 

nonvanishing 03. The residual symmetry is the isospin, 

0(3). 

(ii) The phase with a pion condensate (PIC): the 

phase with the charged pion condensate; |-7T_i_ j 7^ 0. The 
0(3) symmetry is spontaneously broken down to 0(2). 

We now have six GL couplings {a.2,a' 2 ,a.4,, h,0i,02}-> 
and they are functions of thermodynamic variables 
{/1, /ii,T}. We can, in principle, investigate the phase 
structure in this six-dimensional space in full generality. 
However, even if we do that, it would become difficult 
to relate it with the phase structure in the physical pa- 
rameter space. Instead of doing that, here we take the 
advantage of quark loop approximation for which we only 
take into account quark loop effects in the effective poten- 
tial. This would give a reasonable approximation to the 
real potential in the high fugacity region. The feedback 
of quark loops to the potential is 
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where V denotes the spatial volume and Tr should be 
taken over the Dirac, flavor, and functional indices. 
Sq = diag (S u , Sd) is the bare quark propagator, and 
£ = a\ + 123X3 + rfirT 1 is the self-energy for which we 
set 7T2 = without any loss of generality. From this we 
can extract the following explicit expression for a^ in the 
leading order of expansion in jjl\ 16]: 



a4 =ai 0) ( M ,T)+O(^), 
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where we have defined the quantity a 2n for n > 1 as 



^VT).8TA^— — -i 
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with Lu n being the fermionic Matsubara frequency. Per- 
forming the series expansion in jj,i and discarding parts 



of integrand containing a total derivative in p [l6j, we 
can relate 2 and 0i with ai , 



A = Im/^^ + oGu?), 



(4) 



We introduce the GL parameter A by (/3i/mi)Ui=0 
Aa^ ; then we have 
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which we use instead of (/3i//ii)| w= o in the following. In 
order to find explicit expressions for a 2 and a 2 we need 
to specify the model. Let us consider here for a while the 
model with four-fermion interaction of the type 



r- - Q. 



((qq) 2 + (qi~/ 5 Tq) 2 + {qi^qf + (qrq) 2 ) 

+f (fe) 2 + (qiToTq) 2 - {qilaq) 2 - {qrq) 2 ) . 

(6) 
The interaction in the second line comes from so-called 
Kobayashi-Maskawa-'t Hooft determinant, which vio- 
lates the axial U(1)a symmetry explicitly. In this model, 
we have a = — (G + K)(qq), a 3 = — (G — K)(q~T 3 q), and 



a 2 



hc+ck } (l*,T) + O04), 
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TJ ± K +a 2 »>(^T) + 0(rf). 



(7) 



We consider the following two cases here: 

Case (I): Strong U(1) A breaking with K = G. This 
corresponds to the standard Nambu-Jona-Lasinio (NJL) 
model with 



£ int = G((qq) 2 + (qij 5 rq) 2 ). 



(8) 



In this case, a 2 diverges so that one of the chiral four- 
vectors ip becomes irrelevant. In particular, 03 = [3, 

E3- 

Case (II): U(1)a symmetric case with K = 0. In this 
case, «2 = ct' 2 so that we see that the GL potential also 
possesses the symmetry as it should. 

For case I, we find up to the fourth order in fields and 
Mi 

n = ^ + £fV_/ lC7 -4^ 7r i. (9) 

Here, a 2 should be regarded as one evaluated at fii — 
0, i.e., it should be understood as the lowest order in 
expansion in jj,\, Eq. ([7]). If we introduce the notation 



o" + 03 a — a 3 

<Ju = , <Jd : 



2 ' 2 

then we are forced to have a u — ad = a/2 
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FIG. 1. The chiral condensate a (solid), the sigma mass M a 
(dashed, blue), and the pion mass M^ (dotted, red) as a func- 
tion of a2, at /xi — 0. 



For case II, we have at the same order 

n = ^(0 2 + ^ 2 ) 

+4^((^+^ 2 ) 2 +4(0v-(0^) 2 )) (ii) 

(0) ,„n 

-ha + ^^ 2 tt 2 ± + a\ 'X^i(aa 3 + J7V3). 

For stability, we can only investigate the region a 4 > 
at this order, and we have typically A > for \i > 0. 
a u and a a in this case are proportional to uu and dd 
condensates. We notice that the A term in the poten- 
tial contains Aa4/xi(cr 2 — ct 2 ,) so that positive /xi favors 
the situation a ( i > a u . In the following analysis, we sup- 
press the subscript (0) in couplings, so we have /xi and 
four independent parameters {a2, CK4, h, A}, which do not 
depend on /xi. 

Let us first investigate case I. Introducing dimension- 
less parameters a, 7r, fit, and c?2 via 

a = a (V«4) 1/3 , 7r = ^(Va 4 ) 1/3 , 
02 = &2 ft. 2 / 3 a 4 , /ii = /xi (h/ai) 1 / 3 , 

the potential is cast into O = a 4 /i 4 / 3 w, where w does 
not have any explicit dependence on h: 

with (f> = (a, it). When h = and zxj = 0, the system 
has a second-order phase transition at «2 = 0. At finite 
ft ^ 0, the transition gets smoothed to a crossover. In 
Fig.[IJ we show the behavior of chiral condensate a, sigma 
and pion masses (M a , M„) as a function of 0:2- We find 
a crossover from the x^B to an approximately restored 
phase with a ~ when «2 is increased. We can define 
the pseudocritical point 02 = o^ , for example, by the 
point where M a takes the minimum. In this case, it can 
be numerically read as 



o.-, 



1.191 h 2/3 a 



1/3 



At this point, pion and sigma masses are 

M^ c = lMh 1/3 a\ /6 , A/P c = 1.26/i 1/3 a 4 /6 . (13) 

The ratio M% c /Mv c = 1.22 at the pseudocritical point 
is universal to this order of GL expansion. The chiral 
condensate a is also read as 



a = 0.630 ft 1/3 a 4 1/3 (=cr pc ). 



(14) 



For what follows, we concentrate on the effect of /xi at 
the pseudocritical point. In Fig. [3J we show a u , ad and 
\tv± I = 7r as a function of /Xi. We see a second-order phase 
transition to the PIC phase at 



|/*i| = 1-782 /i 1 / 3 ^ 1/3 (= ^). 



(15) 



This can be written in terms of a pc or M% c as 

Ml c = 2.83 a pc = (1.41/^51) M% c . (16) 

We notice that what is in the universal relation in the 
GL framework at this order is the ratio of /xi to the flavor 
singlet quark condensate a rather than that to the pion 
mass M v . Once |/xi| > /Xj, the charged pion condensate 
7T develops. 

A couple of questions are in order here: i) First, one 
might think that /xf determined here by looking at the 
static correlation function might be different from the 
true one, which should be determined by the pole of the 
charged pion propagator due to the kinetic seesaw mech- 
anism [IIIlLSI]- n ) Second, one might wonder what is the 
difference of 7r + and 7r~ condensations. Let us first dis- 
cuss the point i. Actually, this is not the case, and the 
critical point determined by the static effective potential 
exactly coincides with the one by the charged pion prop- 
agators. If we worked out the time-derivative expansion 
in charged pion fields within the Gaussian approxima- 
tion, we would have obtained the effective action density, 
which looks like (with a suitable normalization of 7r fields) 

£ c ff = tr[(i<9 t 7r_L + \pLi,ir±\)^{id t ir±_ + [£i, «"!.])] 

-IM^+TT 2 ), 

where M is the mass parameter being a function of micro- 
scopic/thermodynamic variables, tv± = (wiTi + 7r2T2)/2 
and fi\ = /X1T3/2, with {n, T2, T3} being the Pauli matri- 
ces. The charged pion propagator can be read as 



£^r» = 
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The static part of the propagator is related with the sec- 
ond derivative of effective potential: 
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(19) 
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The determinant of the polarization matrix includes all 
the pole masses: 



(12) 



detDl 



{u 2 - M 2 + ){u 2 - M 2 _) 



(20) 
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FIG. 2. The behavior of condensates as a function of /ii at 
the crossover point, Eq. (|12[) . <r u and <Td(= <r„) are depicted 
by dashed lines (red and blue, respectively), while it is shown 
by a solid line. 



M + = M — Hi corresponds to the tt + pole, while A/_ = 
M + jiti represents the tt~ pole. The critical condition is 
given by the vanishing of either ir + or tt~ mass, and, in 

both cases, M 2 - [i\ = M + M_ = 
This is equivalent with the condition 



det D^j 



= 0. 



d 2 n 



dir 2 



d 2 n 



dir 2 



= 0. 



(21) 



This means that even if the pion masses split due to the 
kinetic seesaw mechanism, the critical chemical point can 
be always obtained by looking at the behavior of static 
susceptibility. Of course, we can only extract from 51 
the multiple of two poles M+M-, not separately M + 
and M-. In other words, it is the geometric average 
of the pole masses that can be read from the curvature 
mass of CI. Let us now come to point ii. Since the time 
dependence of the charged condensate is determined by 



the chemical potential as ir + — ix\ — iit2 



e _Vl *7r, i.e., 

(7Ti , 7T2) ~ (n cos fj.it, 7rsin/xi£), the rotation is counter- 
clockwise when fii > while it is clockwise for /xi < 0. 
We shall refer to the former as 7r + condensation, and the 
latter as 7r~ condensation. 

Now, let us move on to case II. In this case, we have an 
additional GL parameter A > for /i > 0. In the upper- 
panel of Fig. [31 we display the phase diagram in the two- 
dimensional GL parameter space: the (/xi,A) plane. We 
have basically two phases, the xSB and the PIC with a 
7r + or 7r~ condensate. However, in this case, we have 
a u = &d only on two lines specified by \x\ = or A = 
and a u ^ ad in the major part of the plane. This is, 
of course, because of the A term in the GL potential, 
a4\/j,i(a 2 — ct|). On the line A = 0, we have completely 
the same situation as displayed in Fig. [2] in which two 
second-order phase transitions are found at [i\ = ±/ij. 
On the other hand, when A becomes large, the transitions 
eventually change to the first-order ones. We show the 
situation at A = 1 [/i 1 / 3 ^ ] in the lower panel of Fig. EI 
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FIG. 3. (Upper panel): The GL phase diagram in the (/^i,A) 
plane. The solid line represents the first-order phase tran- 
sition, while the dashed line represents second-order one 
separating the x^B and PIC phases. The dashed line in- 
side the shaded region expresses the spinodal line where the 
state without pion condensate becomes unstable. (Lower 
panel): The behavior of a u (dashed, red), ad (dotted-dashed, 
blue), and ty (solid) as a function of /ii at A = 1 h 1 ' 3 ^ . 



We see clearly finite abrupt gaps in the order parameters 
associated with the first-order phase transition. This is 
attributed to the competition between the A term and the 
/ij term in the potential; the former favors a larger \<j u — 
<7d\, while the latter likes the situation n ^ in which a 
smaller \a u — <Jd\ is favorable. Accordingly, there are two 
tricritical points (fii,X) — (±nJ CF , A TCP ) at which the 
second-order phase transitions turn into first order ones. 
Numerically, we find 



m tcp = 2 .02/i 1 /3 a - 1 /3 = 321(Tpc 
A TCP = 0.67/ l 1 /3 a - 1 /3 = 106(Tp ^ 



(22) 



The ratio A TCP / nj cp = 0.33 does not depend on any of 
GL parameters and, thus, is universal to this order. We 
note that the effect of strong flavor mixing due to the 
U(1)a anomaly makes the A term irrelevant via locking 
the condensates a u and ad with the same value [9|, Il7| , 
and. thus, renders the transition a second-order one. 



III. THE GINZBURG-LANDAU APPROACH AT 
SIXTH ORDER 

We now extend the GL analysis up to the sixth order 
so as to explore the influence of /zi near the CP of QCD. 
Doing this in full generality introduces many new GL 
parameters, which makes the analysis quite complicated. 
Instead of doing this, we examine here only case I in the 
previous section, in which one of chiral four-vectors, ip, 
decouples. The general GL for homogeneous condensates 
can be again written in terms of 4> [15| 



n 



-ha + tftf - /3 2 *i 
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rfi 



fti 
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(23) 

GL coefficients a n (n = 2, 4, 6) are expanded in the series 
of fii within the quark loop approximation (up to total 
derivatives) as 



1 a^, 2 b[i\ 
= 10 1 
,0 1 




oci{\x\ = 0) N 
a 4 (//i = 0) 
ae(Mi = 0), 



(24) 



Via explicit computations, we find a = b = 0, and c = 1. 
Similarly, for fiz, /Ji's we find 
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{/34,/3 46 } 



= M? 
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(25) 



with d = -1/4, e = 0, and / = f b = -2. Putting 
them all together, and assuming the condensate to have 
a charged pion component, i.e., <fr = (a, n, 0, 0), we arrive 
at 



f2 = —ha 



&2 „2 



2 ~ 



Oi+JJ^ae 4 1 04-Hf£«6/2 2_2 . a 4 4 
4 a ^ 2 a n + ~T W 
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Here «2j <*a an d ocq should be understood as those eval- 
uated at [i\ — and, thus, are functions of fj, and T only. 
a e has a dimension of energy so we use this as the 
unit of energy. In the following analysis, we, thus, set 
a& = 1, but the proper dimension of any quantity can 

— 1/2 

be recovered any time by use of a e . Now, assuming 
h > 0, we get rid of h via scaling 



a = ah 1 / 5 , ir^irh 1 / 5 , W =/i I /i 1 / 5 , 



a-2 



a 2 



h"l\ 



o?4 = 0.4 



h 2 'K 



(27) 



h & ' b uj. 



Then, h is scaled out of the potential as Q. 
Therefore, we now need to explore the GL phase diagram 
in the three-dimensional GL parameter space (a<2, a.4, /ii). 
How does (0:2,04) map onto the (fJ.,T) -plane? Before 
going into the discussion of GL phase diagrams, let us 
briefly sketch how the GL parameters (0:2, 0:4) map onto 
the (n, T) plane taking the chiral limit (h — 0) for sim- 
plicity. Figure S] illustrates how (02,04) spans the local 
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FIG. 4. The illustrative figure that shows how (0:2,0:4) spans 
the local coordinate in the vicinity of the TCP within the NJL 
model defined by Eq. |(HJ and the three-momentum cutoff A. 
This is depicted for two values of coupling g = GA 2 = 2.0 
(upper panel) and g — 2.0 (lower panel). The solid line (red 
online) shows the curve on which 02 = 0, while the dotted- 
dashed line (blue online) does that for 04 = 0. The point 
of intersection gives the location of the TCP. The region for 
04 > is shaded. The solid line in the shaded area represents 
the second-order chiral transition, while that in the unshaded 
area only gives the spinodal line on which the Wigner phase 
with £7 = ceases to be a local minimum. 



coordinate in the NJL model specified by Eq. © and 
three-momentum cutoff A. The upper panel shows the 
situation for g = GA 2 =2.5, while the lower panel shows 
that for g = 2.0. In the figure, the solid line expresses the 
curve determined by 02 = 0, which separates the (n,T)- 
plane into two regions, one for 02 > and the other for 
o 2 < 0. Similarly, the dotted-dashed line shows the curve 
for 04 = 0. The point of intersection determines the lo- 
cation of the TCP. The region 04 > is shaded just for 
a guide. The solid line inside this region determines the 
second-order chiral phase transition, while that outside 
it only specifies the spinodal line on which the Wigner 
phase ceases to be even a local minimum. The axises 
of the local coordinate system (02,04) are depicted by 
arrows starting from the TCP; the 02 (04) coordinate 
points to the direction for 02 > (04 > 0) side with 
being tangent to the line of 04 = (02 = 0). The loca- 
tion of the TCP as well as how (02,04) maps onto the 
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FIG. 5. The GL phase diagrams in the absence of in. The 
phase diagram in the chiral limit h — (upper panel) and 
the one off the chiral limit (lower panel). The shaded area 
represents the spinodal region. See the text for detail. 



(/i, T) plane depends on the detail of the model or reg- 
ularization scheme, so, for what follows, we only discuss 
the phase diagram in the (ev2,ev4) plane. 

GL phase diagram for quark matter without an isospin 
imbalance: Let us start with the case with /ii = h = 0. 
In the upper panel of Fig. [5j the phase diagram for this 
case is displayed. In this case, h in the figure labels can be 
regarded as an arbitrary unit; changing it does not mod- 
ify the phase diagram. There are two phases: the xSB 
phase with a ^ and the symmetric (Wigner) phase. 
For CV4 > 0, these phases are separated by a second-order 
phase transition located at a 2 — 0, which is depicted by 
the dashed line. For en < 0, the transition is replaced 
by the first-order one at c*2 = j§oi\ shown by a solid 
line [lg|. Accordingly, the TCP is located at the origin. 
The shaded area shows a spinodal region in which one 
of the states exists as a local minimum of the potential; 
at eti — for a a < 0, the Wigner phase becomes un- 
stable, while at the line a-i = \oi\ for ct± < 0, a local 
minimum corresponding to the x>3B state vanishes. In 
the lower panel of Fig. [5j we show the phase diagram for 
h y^ 0. In this case, we have only the situation a ^ 0. 
Nevertheless, the first-order phase transition survives and 



separates the xSB phase and a nearly symmetric phase 
with cr^O. The first-order phase transition ends at the 
CP. The exact location of the CP, (a§ F t&i): is derived 
analytically [2(| 



/„CP „CP\_/5 3 4/5 l4/5 5 1.2/5^ 

("2 i a A ) — ^42275" ' _ 2 1 /533/5 ft J 



(28) 



This is numerically evaluated as ~ (2.28/i 4/5 , ™2.25/i 2/5 ) 
The chiral condensate at this point is found as 



= 0.822 h 1/5 (=<j ). 



(29) 



The dotted line starting from the CP expresses the pseud- 
ocritical line determined by the location of the minimum 
in the sigma meson mass. The shaded area again repre- 
sents the spinodal region, in which there is another state 
competing with the ground state. 

GL phase diagram for quark matter with an isospin 
imbalance: Now, we discuss the influence of \i\ on the 
phase diagram. The phase diagrams for several finite val- 
ues of Hi are shown in Fig. [5] From Fig.^a) to Fig.[BJd), 
the value of /zt increases. We now suppress the spinodal 
region surrounding the line of the first-order chiral phase 
transition. 

Let us start with (a), where the value of [i\ is set to 
0.2ft 2 / 5 , that is, in terms of Co, /ii = 0.54ero- In this 
case we notice that the structure in the proximity of the 
CP is unaffected even though the location of the CP is 
shifted downward according to 



/<r 



(a«>4 CP )^(«?>? P -Mi 2 )- 



(30) 



This can be easily understood from the coefficient of the 
cr 4 term in Eq. (|2l)| . Recalling the direction in which the 
local coordinate a 4 points in the (^,T) plane (Fig. 01 , 
we expect that the CP shifts to the direction of lower 
temperature and higher chemical potential due to the 
isospin asymmetry. This is consistent with the analysis 
done in Ref. [21(. Moreover, a recent analysis within a 
specific model shows at some critical value of [i\ the TCP 
(CP) can even touch the 11 axis disappearing eventually 
from the phase diagram [22]. Apart from this simple shift 
of chiral phase transition, we notice that the continent of 
the PIC dominates the region deep in the xSB. The PIC 
and xSB phases are separated by a second-order phase 
transition in all regions shown in the figure. 

Figure [H^b) shows how the situation changes when H\ 
is increased to /i 2 = 0.5/i 2//5 , that is, fii = 0.8600. We 
notice that the CP moves further downward according 
to Eq. (|30p . and the continent of the PIC gets wider as 
expected. Moreover, the transition from the xSB to the 
PIC now has a branch of the first-order phase transi- 
tion, which is drawn by a solid line. Accordingly, there 
appears a new tricritical point denoted by TCP' on the 
critical line. This is actually the tricritical point at which 
three critical lines meet up, once we introduce an exter- 
nal field for the charged pion condensate. The line for the 
first-order phase transition departing from the TCP' en- 
counters the line of the first-order chiral phase transition 
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FIG. 6. The GL phase diagrams for several values of /xi; from (a) to (d), fii increases as (a) p 2 — 0.2/i 2 ' 5 , (b) <x 2 — 0.5/i 2 , 
(c) nf = (fif) 2 , and (d) /j, 2 = 3.0/i 2/5 . For nf, see the text. 



at the point "TP" , which stands for the triple point. At 
the triple point, three phases — the xSB, a nearly restored 
phase, and the PIC — coexist and compete. 

Now, let us discuss Fig. |6]Jd) before Fig. [6jc). In 
Fig. E^d), Hi is set to a large value ttj = 3.0/i 2 / 5 , which 
corresponds to i*i = 2.1<r - hi this case, the PIC domi- 
nates a major part of the phase diagram, and the struc- 
ture of the chiral phase transition is now completely hid- 
den. There remains a dotted line outside the PIC, which 
is just the chiral crossover. The transition from the PIC 
to the xSB phase becomes widely of first order, and, ac- 
cordingly, the location of the TCP' is shifted upward. 

Since the CP is completely hidden in Fig. [SJd), there 
should be a critical value of /ii at which the CP vanishes 
from the phase diagram. In Fig. Hfc), we show this situ- 
ation. The critical chemical potential is /ij = lA77h 2 ^ 5 , 
which translates into /ii = 1.48<to = Mi- As shown in the 
figure, the CP comes across the line expressing the first- 
order phase transition between the PIC and %SB phases. 
This means when fii becomes large, the chiral critical 
point could fade out from the QCD phase diagram. The 
ratio of critical /xi to &$, the chiral condensate at the CP, 
can be numerically evaluated as ~ 1.48, which is uni- 
versal being independent of any GL parameters to this 
order. 



IV. CONCLUSION 

We performed a systematic GL analysis on the effect 
of isospin asymmetry on the chiral crossover, the CP and 
its neighborhood. We first focused on how the crossover 
is affected by the isospin density. To this aim, we de- 
rived a general GL potential up to the quartic order in 
two chiral four-vectors, 4> and its parity partner </?. Mak- 
ing use of the quark loop approximation together with a 
perturbative expansion in fj,j, we have studied not only 
the nature of the phase transition to the PIC but also 
how it is affected by the effect of the U(1)a anomaly. We 
found the effect of the isospin-flip odd A term in the po- 
tential makes the phase transition to first order at large 
A. Since A vanishes at fi = and increases with /u, this 
may explain why the transition to the PIC is observed to 
be first order at finite /x in several model analyses [y, [7| ■ 
The effect of flavor mixing due to the U(1)a anomaly 
was found to diminish the effect of the A term by lock- 
ing two condensates, a u and ad- We have derived three 
model-independent universal ratios — /iJ/(Tp C at A = 0, 
A p /<7 pc and /m[ cp /a pc at the TCP — , which are inde- 
pendent of any GL parameters to the fourth order. 

We then extended the analysis up to the sixth order 
of GL expansion, so as to study the isospin effect on 
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the CP. Restricting the analysis to the case with strong 
U(1)a symmetry breaking, we studied how the CP and its 
neighborhood are affected by the incorporation of isospin 
density. We found that it has remarkable effects; it not 
only causes a shift of the location of the CP, but also 
brings about the development of a sizable region for the 
homogeneous pion condensate. This leads to the appear- 
ance of new tricritical and triple points. Moreover we 
showed that the CP disappears once the isospin chemical 
potential is increased above a critical value. We derived 
the critical value /Xj and a universal relation relating it 
with the size of chiral condensate at the CP, <To . 

There are several directions into which the current 
work can be extended. First, we need to take into ac- 
count the possibility of inhomogeneous phases since they 
are known to play an important role near the TCP/CP 
(m 1231-2 8]. This is actually now under investigation [14j . 
Second, the extension to three flavors would be inter- 



esting. This would require the incorporation of a kaon 
condensate and a diquark condensate of the color-flavor 
locked type. In particular, it is known that the inter- 
play between the chiral and diquark condensates via the 
axial anomaly leads to a rich variety of phases and an 
appearance of new multicritical points [23 . |30|. Third, 
the effect of vector interaction should be taken into ac- 
count [3l|, [32| . This is especially needed when we look 
at dynamical aspects of the critical behavior. In fact 
the dynamic universality class of the chiral CP is known 
to be the same as the liquid-gas CP [33|, [3J] . Lastly, it 
is strongly desirable to seek the stiff experimental sig- 
natures of critical points observed here, such as those 
discussed for CP [35j . 

The author thanks K. Suzuki for several useful com- 
ments. A part of numerical calculations was carried out 
on SRI 6000 at YITP in Kyoto University. 
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